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Q . Abstract. In order to include a correction by the Coulomb interaction in Bose-Einstein correlations 

' (BEC), the wave function for the Coulomb scattering were introduced in the quantum optical 

approach to BEC in the previous work. If we formulate the amplitude written by Coulomb wave 
■ functions according to the diagram for BEC in the plane wave formulation, the formula for Stt^BEC 

becomes simpler than that of our previous work. We re-analyze the raw data of Stt^BEC by NA44 
" ^ . and STAR Collaborations by this formula. Results are compared with the previous ones. 
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; INTRODUCTION 

O ; 

^ . Recently, in addition to tlie data on the charged two-body Bose-Einstein correla- 

PhI tionAiBEC), data on the three-body charged BEC have been reported IB S S 0]- hi 

p • some papers 10, 111], Coulomb correction is done with fixed interaction region, 5 fm. 

1^ . On the other hand, raw data with acceptance correction are reported lIlL0]- In Ref. fs'], 

\ the 3;r BEC is calculated by the use of the Coulomb wave function with fixed source 

\ radius R. We have analyzed the BEC without assuming the size of source radius using 

^ ■ the Coulomb wave function |@,0]- In this paper, we would like to refine the formula 

^ ! for the 3;r BEC using the diagrammatic representation for BEC in the quantum optical 

(QO) approach [8, 9, 10]. 

In order to describe the two-body charged BEC ( for example, 2;r^ system ), we should 
solve the Shrodinger equation of Coulomb scattering. The regular solution at the origin 
of the Coulomb potential is given by, 

V^k,(x,,x,) = r(l +/77,,0e-^'''v/2e*o-roF[-/7],,, 1; (fc„r,-, -k,-, -r,-,)], (1) 

where, r,-^- = x, - xj, k/j = (mjk,- - m/kj)/ (m,- + m^), rtj = |r/j|, kij = \kij\ and riij = 
eie jUij/kij. mjj is the reduced mass of m,- and nij, i,j = 1,2,3 and i ^ j. F[a^b\x] is the 
confluent hypergeometric function, and T{x) is the Gamma function. 



The wave function of identical bose particles should be symmetrized. For two particle 
momentum density we have, 



N^^'''^ = -n / P(x/)^^x,|v/^^k2(xi, X2) + <k2(x2, xi; 

2 

= Yl p{x,)d\,{Gi+G2), 

Gl = ^(|v^^,k2(Xl, X2)|^+|V^^jk2(X2, xi) 

G2 = Re(v/^^k2(xi, X2)v/fk2(x2, xi)), 



where, 



The three particle momentum density for 37t~ BEC is written as 

^^'"'^ = ^tU [ Pi^.)d'x,\l^A{j)\\ (3) 

A(l) = Ai = V/^jk2('^l' X2)v4k3('^2, X3)i/4ki(x3, Xi), 

A(2) = A23 = V^,k2(xi' X3)v4k3('^3, X2)V/^3ki(x2, Xi), 

A(3) = A12 = V^k,k2(x2, xi)v4k3('^i' X3)v/^3kj(x3, X2), 

A(4) = A123 = V^ik2('^2, X3)v4k3(x3, xi)v/£k,(xi, X2), 

A(5) = Ai32 = V^kik2N, xi)v4k3('^i' '^2)V^,k,(x2, X3), 

A(6) = Ai3 = v/^jk2(x3, X2)v4k3('^2, xi)v/^3k;(xi, X3). (4) 
In the plane wave approximation, each amplitude A (/) approaches to the following form; 

A(l) = Ai e'l'i2-"'i2g!k23-r23g''k3rr3i ^ g(3/2)/(ki-xi+k2-X2+k3-X3) 

A(2) = A23 e*'2'''i3e''*23-r32g'k3rr2i ^ g(3/2)/(ki-xi+k2-X3+k3-X2)^ 

^(•3^ _ g'ki2-r2ig;k23-ri3g;k3rr32 _ g(3/2);(ki-X2+k2-xi+k3-X3) 

A(4) = A123 e''*i2-"'23g'k23-r3igik3i-ri2 ^ g(3/2)/(ki-X2+k2-X3+k3-xi)^ 

^(^5-) ^ ^^22 Z!^ g'ki2-r3ig!k23-ri2g;k3i-r23 ^ g(3/2);(ki-X3+k2-xi+k3-X2)^ 

A(6) = Ai3 e*'2-r32g!k23-r2i^/k3i-ri3 ^ ^(3/2)/(ki-X3+k2-X2+k3-xi)^ ^^-^ 

where PW means the plane wave approximation of the amplitude. It should be noted 
that factor 3/2 does not appear in the formulation of 3;r^BEC using the plane wave. 



QUANTUM OPTICAL APPROACH 



The amplitude squared in Eq. @ can be classified into the following groups, 

Fi = (l/6)[AiAl+Ai2Al2+A23A53+Ai3Al3+Ai23Al23+Al32Al32], 

Fn = (l/6)[AiA|2+A23Al23+Al3Al32 + C.C.], 

F23 = (l/6)[AiA^3+Ai2At32+Ai3At23 + C.C.], 

F31 = (l/6)[AiAt3+A23At32+Ai2At23 + C.C.], 

^123 = (l/6)[AiAt32+Ai32At23+Ai3At2+Ai2A^3+A23At3+Ai23At] 

i^l32 = (l/6)[AiAt23+A23At2+Ai2At3+Ai23At32+Ai32At+Ai3A^3], (6) 

where, c.c. denotes the complex conjugate. In the plane wave approximation, F\ reduces 
to 1, Fij corresponds to exchange between i and j charged particles, and F123 correspond 
to exchange among three charged particles. 

In the previous work [6], we introduced the Coulomb wave function to the core-halo 
model [Tlj, where following parameters are used; the fraction of multiplicity from the 
core part, fc = {ricore) I (ptot), and the fraction of coherently produced particles from the 
core part, pc = {uco) / {ncore)- In the quantum optical approach, chaoticity parameter is 
defined as ;? = {richao) / {f^core) , where p = I — Pc- In the following, we use fc and p. The 
radius of halo part is assumed to be infinitely large. Therefore, particles emitted from the 
halo part do not contribute to BEC, namely ex(^ange of particles. In Ref. LIOL the higher 
order BEC is formulated, where a contamination effect is included. If we can identify 
the contribution of the halo part to the contamination, both formulations coincide. 
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FIGURE 1. Cumulant up to the third order, (/,;) = (2,3) or (3,2). 

In the QO approach, the two particle momentum density p(pi,P2) and the three 
particle momentum density p(pi, p2, P3) are given respectively as [.9.. .10,1. 

P(Pl,P2) = P(Pl)p(P2)+^(Pl,P2), 
P(Pl,P2,P3) = P(Pl)p(P2)p(P3)+^(Pl,P2)p(P3) 

+,?(P2,P3)p(Pl) +.?(P3,Pl)p(P2) +^(Pl,P2,P3)- 



The second order cumulant ^(pi , P2) and the third order cumulant g(pi , P25 Ps) are given 
in Fig.[T][9, 10]. In Fig.Ul the chaotic component r(pi,p2) is shown by the solid line 
with arrow from 1 to 2, and the coherent component c(pi,p2) is shown by dotted line. 
In the formulation of BEC with the Coulomb wave function, momentum and coordinate 
cannot be decoupled. Therefore, interpretation of diagram is somewhat modified; The 
source function p (x) is attached to the starting point of solid line, and the delta function 
d{x) is to that of dashed line. In addition, p is replaced by fcp, and I— p ishy fc{l— p) 
to include the contribution from the halo part. 

In the 2;r"BEC, N^^ corresponds to p(pi,p2), Gi to p(pi)p(p2), and G2 to 
g(pi, P2) beside the normalization factor. Then, the formula for 2;r^BEC is written as, 

2 . 

Cll pixi)d\{Gi+f^p^G2) 

+C J d\i J j3x2p(xi)5(x2)2/2p(l-p)G2 

In the third order BEC, A^-^'^ corresponds to p(pi,p2,P3), Fi to p(pi)p(p2)p(p3), 
Fij to ^(p/,P;), and (F123 +F123) to ^(pi,P2,P3)» beside the normalization factor. Then, 
the formula for 3;r^BEC is given by 

]^ = Cllj p{x,)d\i[F,+3f^p^Fn + 2f^p'-Re[Fn3]] 

3 f 

+ Cll d\,5{xi)p{x2)p{x3)6f^p{l-p){Fn+fcpRe[Fn3]). (7) 
Compare Eq.lQ with Eq.(12) in Ref. where F2 = F12 + F23 +F31 and F3 = F123 + 

-P'l32- 



ANALYSIS OF 3 tt BEC 

The formula Q is applied to the analysis of raw data on 3;r BEC by STAR Collabora- 
tion [,4J and NA44 Collaboration fl]. 



TABLE 1. Analyses of raw data on 3;r" BEC by STAR Collaboration jEJ]. 





P 


1.0 


0.8 


0.6 


0.40±0.02 


present work 


fc 


0.75±0.02 


0.78±0.02 


0.85±0.02 


l.O(fixed) 




R{fm) 


5.34±0.24 


6.21±0.30 


6.88±0.35 


7.43±0.34 






2.68/34 


0.93/34 


0.65/34 


0.58/34 


Ref. \6} 


fc 


0.75±0.02 


0.78±0.02 


0.85±0.02 




R{fm) 


5.34±0.24 


5.99±0.28 


6.48±0.32 








2.80/34 


1.39/34 


0.98/34 





The results for 3;r BEC by STAR Collaboration [Q] are shown in Table[T]and Fig.|2l 
For comparison, the results in the previous work are also shown in the lower part of 
Table [B Parameters fc and p are almost the same with the previous calculations. The 



radius R is almost the same at p = 1.0, and the present result gradually becomes larger 
than the previous one. At p = 0.6, the new result becomes about 6% larger than the 
previous one . Fitting to the STAR data is slightly improved. In Fig|21 the result at p=l 
is shown, and the possible region of p and fc estimated from the analysis of 2;r^BEC 
and 3;r^BEC is shown. Each hatched region is defined by the fitted value ± 2 standard 
deviation. 



1.8 
o 1.6 
5 1.4 
^1.2 
1 

0.8 



' ' ' ' I ' ' ' ' I ' ' ' ' I 
STAR Au+Au 130 A GeV" 
p = 0.40 + 0.02 
^ = 7.4 + 0.3fm J 




50 100 150 
Qi (MeV) 



200 




FIGURE 2. (a) Analysis of raw data on 3;r"BEC by STAR Collaboration; is fixed at 1 . (b) Possible 
region of fc and p estimated from the analyses of In^ BEC and 37r^BEC. 

TABLE 2. Analyses of raw data on 37t- BEC by NA44 Collaboration f]^. 





P 


1.0 


0.8 


0.6 


0.32±0.05 


present work 


fc 


0.67±0.04 


0.70±0.04 


0.76±0.05 


l.O(fixed) 




R{fm) 


2.88±0.39 


3.32±0.49 


3.66±0.54 


4.05±0.57 






6.64/15 


6.78/15 


6.81/15 


6.80/15 


Ref. [TJ 


fc 


0.67±0.04 


0.70±0.04 


0.76±0.05 




R{fm) 


2.89±0.39 


3.22±0.46 


3.47 ±0.51 








6.7/15 


6.8/15 


6.8/15 





o 1 

pa 
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FIGURE 3. (a) Analysis of 37r-BEC by NA44 Collaboration; is fixed at 1 . (b) Possible region of fc 
and p estimated from the analyses of 2n^ BEC and 37r^BEC. 

The results for 3;r^BEC by NA44 Collaboration 1 1] are shown in Tableland Fig.|3] 
Previous results [7] are also shown in Table |2l Parameters fc and p are almost the same 
with the previous calculations. As for the radius R, the present result gradually becomes 
larger than the previous one, as the chaoticity parameter decreases from 1 to 0.6. At 



p = 0.6, the new result becomes about 5% larger than the previous one. Fitting of the 
present calculations to the NA44 data is almost the same with that of the previous results. 



CONCLUDING REMARKS 

We refine the formula for 3;r BEC according to the diagram in the QO approach. In 
our formulation three parameters, fraction of core part chaoticity parameter p, and 
radius of interaction range R are included. We apply the formula to the analysis of 
data on 3;r^BEC by STAR and NA44 Collaborations, and compare with the previous 
calculations. Parameters fc and p are almost the same with the previous calculations. 
The radius R is almost the same at p = 1.0, and the present result gradually becomes 
larger than the previous one, as p decreases. At p = 0.6, the new result becomes about 
5% larger than the previous one. Fitting to the STAR data is slightly improved. That to 
the NA44 data is almost the same. 

From the comparison of the results on 2;r~BEC with those on 3;r~BEC, we 
have a relation. 



^2nd 



l.5R3rd- (8) 



Factor 1 .5 in Eq. dH) corresponds to 3/2 in Eq. (O. This result suggests that our approach 
should be equivalent to the plane wave formulation applied to the data with Coulomb 
correction. Detail calculation is reported in Ref. Il2l] 
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